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Abstract 

The effect of acoustic coupling on random and harmonic plate vibrations is stud- 
ied using two numerical models. In the coupled model, the plate response is obtained by 
integration of the nonlinear plate equation coupled with the nonlinear Euler equations for 
the surrounding acoustic fluid. In the uncoupled model, the nonlinear plate equation with 
an equivalent linear viscous damping term is integrated to obtain the response of the plate 
subject to the same excitation field. For a low-level, narrow-band excitation, the two models 
predict the same plate response spectra. As the excitation level is increased, the response 
power spectrum predicted by the uncoupled model becomes broader and more shifted to- 
wards the high frequencies than that obtained by the coupled model. In addition, the 
difference in response between the coupled and uncoupled models at high frequencies be- 
comes larger. When a high intensity harmonic excitation is used, causing a nonlinear plate 
response, both models predict the same frequency content of the response. However, the 
level of the harmonics and subharmonics are higher for the uncoupled model. Comparisons 
to earlier experimented and numerical results show that acoustic coupling has a significant 
effect on the plate response at high excitation levels. Its absence in previous models may 
explain the discrepancy between predicted and measured responses. 



1. Introduction 

It is well known in structural dynamics that linear plate theory cannot be used 
to accurately predict structural responses at high excitation levels. In an attempt to over- 
come this weakness, several nonlinear plate models have been proposed. The nonlinearities 
introduced in the various models can be classified in two categories: geometric nonlinearity 
or material nonlinearity. The nonlinear model based on the geometric nonlinearity, also 
referred to as the large deflection model, has been extensively used in the literature. 1-10 
Although this model gives better predictions than those of the linear theory, it overestimates 
the frequency content or “broadening” of the response spectrum at high excitation levels. 

To further improve the predictions of the nonlinear analysis, several damping mod- 
els have been studied. Linear, nonlinear, and viscous damping are among the models 
used. 11-15 Prasad and Mei 16,17 used nonlinear damping in their large deflection model. 
A set of modal equations of the Duffing form coupled in the nonlinear stiffness and uncou- 
pled in the nonlinear damping was derived. An approximate solution was obtained for rms 
quantities (such as displacement and strain) and spectral density functions by the equivalent 
linearization method. The nonlinear damping was found to contribute to the broadening 
of the power spectra at high excitation levels. However, Moyer 1 ® and Reinhall and Miles 10 
showed that the method of equivalent linearization gives inaccurate response spectra for 
large deflections. The broadening of the response spectrum at high excitation levels was 
found to be caused by the nonlinear stiffness. 19 Robinson and Mei 20 studied the influence of 
nonlinear damping on panels random response using a time domain simulation. They found 
that the nonlinear damping was responsible for a narrowing of the nonlinear response power 
spectrum. 

Recently, Robinson etal. 21 used a finite element method to integrate the nonlin- 
ear plate equations with an equivalent linear viscous damping term to study the linear and 
nonlinear response of a plate subjected to a narrow band random excitation. Their results 
were in good agreement with experiments for low excitation levels. However, for high excita- 
tion levels, the model overpredicted the broadening of the response spectrum and predicted 
higher response levels at high frequencies. 

In most structured dynamics analyses, there is little discussion about the effect 
of the surrounding acoustic fluid on the response of a structure to a given excitation field. 
Prendi, etal. 22,23 showed that a strong coupling between plate vibration and the surrounding 
flow fields exists at high excitation levels. 

In this paper, an attempt is made to explain the discrepancies between the experi- 
mented and numerical results obtained by Robinson, etal. 21 Using the model developed by 
Frendi, etal. 23 and changing the configuration of the computational domain to match that 
of the experiments, the response of a flexible plate to both a narrow band random excitation 
and a harmonic excitation is studied using two models. One model accounts for the full 
nonlinear coupling to the surrounding fluid. The other replaces the nonlinear coupling by 
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an equivalent linear viscous damping on the plate. 

The remainder of the paper is organized as follows. In section 2, a detailed descrip- 
tion of the analytical model is given. Section 3 describes the numerical techniques, and the 
results are discussed in section 4. Finally, the conclusions are given in section 5. 


2. Formulation of the model 

As shown on fig. 1, the computational domain is composed of three regions. The 
acoustic fluid in the cavity region (or Top domain) is separated from that in the open space 
region (or Bottom domain) by a flexible plate clamped between two rigid plates. Acoustic 
disturbances are introduced at the left boundary of the Top domain and propagate over 
the rigid and flexible surfaces. There is no meanflow on either side of the flexible plate. 
The governing equations in the acoustic fluid regions are the two-dimensional, compressible, 
nonlinear Euler equations. In a cartesian coordinate system, x and y, these equations can 
be written in conservation form as 


dQ_ = dF_ dG 
dt dx dy ’ 


( 1 ) 


where Q is the vector ( p , pi t, pv, e)^, p is the density, pu and pv are the x and y momenta 
respectively, and e is the toteil energy per unit volume given by 


e = ^p(u 2 + v 2 ) + P<h)T. 
In eq. (1), the functions F and G are: 

( P u > 

F= pu 2 +p 

puv 1 

\«(e + p)/ 

and 

f pv y 

puv 

vj — 2 . 

PV +P 
\u(e + p) ) 

In addition to Eq. (1), an ideal gas state equation is used 


( 2 ) 


(3) 


P = PRT, 


(4) 


where p is the pressure, p the density, R the gas constant, and T the temperature. 

The equation describing the motion of the one-dimensional flexible surface is 

d^w d 2 w d 2 w . dw 

D ~=~ r ~ Nx ~fa2 + p P h ~^2 + + r °)'S7 = ^ Pt 


dx* 


at 2 


dt 


(5) 
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where w is the plate transverse deflection, p p the mass per unit volume of the plate, and h 
the plate thickness. The total viscous damping is the sum of the structured damping, T,, and 
the acoustic damping, r a . When the plate vibration is coupled to the surrounding acoustic 
fluid, T a = 0. In eq. (5), D — Eh 3 / 12(1 — v 2 ) is the stiffness of the plate, with E being the 
modulus of elasticity and v the Poisson ratio of the plate material. The coefficient AT* in eq. 
(5) is given by 

Eh f x o+ L (dw\ 2 

~ 2 L Jx o \dx) 


N x 


dx, 


( 6 ) 


which represents the tension created by stretching of the plate due to bending. In eq. (6) x$ 
is the origin of the flexible plate and L its length. The forcing term on the right-hand-side 
of eq. (5) is 

Ap = p~-p + , (7) 


where p* and p~ are the pressures on the surfaces above and below the plate, respectively. 
When acoustic coupling is neglected, the plate response is obtained by neglecting p in eq. 
(5) and using the same p + at the plate center as that obtained by integration of the Euler 
equations. This is done in order to reproduce the results obtained by Robinson, eial, 21 

Equations (l)-(7) are written in a nondimensional form using the following reference 
quantities for the different variables: 


(s,y»™)re/ = l ref> 


tref — 



) 

Cqo 


and 



Pref = Poo, ( u >v,-^-) re f = c oo, and (p, e) re f = Poocjp. (8) 

The notation (•, -) re f is used to represent the reference quantity of the variables in parenthe- 
ses. The various freestream fluid properties are those of air at sea level conditions which are: 
temperature Too = 519 °R, density poo = 1.147xl0 -7 lb sec 2 /in 4 , pressure poo = 14.7 psi and 
sound speed Coo = 13392 in/sec. The specific heat at constant volume is Cy = 6.1776xl0 5 
in 2 /(sec 2 °R), the ratio of specific heats is 7 = Cp/c* = 1.4, and the reference length is 
l re f = 12.0 in. 


3. Method of solution 

The unsteady Euler equations (eq. (1)) are solved using an explicit finite difference 
scheme. The scheme, which is a generalization of MacCormack’s scheme obtained by Gottlieb 
and Turkel, 24 is fourth order accurate in space and second order accurate in time. The 
numerical scheme, applied to a one-dimensional equation of the form 


du dF 
dt dx ’ 

consists of a predictor step given by 

u*i=u? + ^(-7 Fi + 8 F i+1 - F i+2 ), 


( 9 ) 


( 10 ) 
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followed by a corrector step of the form 


n-j-1 __ £ 


Ui = 


2 L 


u*i + u? + ^(7 Ft - SFtt + F?_ 2 ) 


( 11 ) 


In the above equations, the subscript i denotes the spatial grid point and the superscript n the 
time level. The fourth-order accuracy is obtained by alternating the scheme given above with 
its symmetric variant. 2 ^ Operator splitting is used to reduce the two-dimensional problem 
to a sequence of one- dimensional problems. If L x and Ly denote the solution operators for 
the one-dimensional x and y problems, then the solution to eq. (1) is obtained by 


Q n+2 = L x LyL y L x Q n . 


( 12 ) 


Further details regarding the method and the advantage of fourth order schemes can be 
found in Bayliss, etal. 25 

The boundary conditions employed on the rigid and flexible surfaces for the Euler 


equations are 



v = 0 

and T — 

(13) 

over the rigid surfaces, and 



dw 

V= ~di 

and T = T W 

(14) 


over the flexible surface. In eqs. (13) and (14) T w is a specified wall temperature; in this 
paper it is taken to be the free-stream temperature, Too- The x-component of the velocity 
(u) is obtained through linear extrapolation from the interior over both the rigid and flexible 
surfaces. 

The pressure boundary conditions are as follows: over the rigid surfaces where the 
time rate of change of the normal momentum is zero, the pressure is calculated using the 
normal momentum equation by simply imposing the normal gradient of the sum of pressure 
and vertical momentum flux to be zero (gy(p + pv 2 ) = 0). Over the flexible part of the 
surface, a linear extrapolation from the interior is used. At the left boundary of the Top 
domain, a nondimensional perturbation velocity is specified ti as 


u = eR(t) or u = esin(u;t), 


(15) 


where R(t) represents a narrow-band random excitation and u is a nondimensional pertur- 
bation velocity. In eq. (15), e is the peak amplitude of the disturbance for the harmonic case 
(with u) being the frequency), and an rms amplitude in the random case. Using the routine 
RNNOF of the IMSL library 26 , Gaussian random data is generated; then, using a quadratic 
narrow-band filter, the data is filtered to give the desired input. The pressure and vertical 
velocity (v) are obtained by linear interpolation from the interior, and the temperature is 
imposed to be constant (T = T w ). The remaining non-physical boundary conditions (Top 
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domain right boundary, Bottom domain right, left, and bottom boundaries) are derived using 

« • 07 

the method of characteristics. 

The plate equation is integrated using an implicit finite difference method for struc- 
tured dynamics developed by Hoff and Pahl. 2 ® The calculation of N x weis done using Simp- 
son’s rule of integration. The boundary conditions used to solve the plate equation are those 
for a clamped plate 

dxv 

w — -q— = 0 at x — xq, xq + L. (21) 


The coupling between the acoustic fluid and the flexible plate is performed as 
follows. The method used in the fluid is explicit. Thus the pressure fields in the Top 
and Bottom domains are updated by using the value of ^ at the previous time-step as 
a boundary condition for the Euler equations. Then, using the new values of p"*" and p~ , 
the plate equation is integrated to obtain the new vertical velocity This procedure is 
repeated at every time-step. 

For every coupled calculation, the time history of the pressure at the center of the 
flexible plate in the Top domain is stored in order to be used as input for the equivalent 
uncoupled case. This was done in order to match the numerical procedure used by Robinson, 
eta l. 21 


4. Results and Discussion 

The results presented in this paper are obtained for a flexible plate having the 
following properties: stiffness D = 1095.6 lbf • in, mass per unit area p p h = 2.21xl0 -5 lbf 
• sec 2 /in 3 , and a Poisson ratio v = 0.3. Two values of total viscous damping (r a + T 0 ) 
are used; 2.5xl0 -3 lbf • sec/in 3 used in the fully coupled calculation (where T a — 0) and 
5xl0 -4 lbf • sec/in 3 used in the uncoupled calculation. These two values axe chosen such that 
the two calculations predict the same response for low-intensity excitations (linear vibration 
regime). The plate is 15 in long, 11 in wide, and 0.13 in thick. The first natural frequency 
of the plate is 112 Hz. The Top domain of Fig. 1 is 12 in high in the y-direction and 120 
in long in the x-direction, and the number of computational points used are 101 and 181 in 
both directions, respectively. The Bottom domain is 120 in long in both directions with 161 
points used in each direction. 

The power spectra shown here are for the center plate displacement response and 
the pressure on either side of the plate center. Also typical instantaneous pressure distribu- 
tions are shown for both the Top and Bottom domains. 

• Previous Experimental and Numerical Results 

At first the results obtained in Ref. 21 are summarized. Figure 2 shows the power 
spectral density of the strain obtained both experimentally and numerically for a narrow- 
band excitation sound pressure level of 140 dB. The figure shows that the power spectrum 
obtained from the simulation is in good agreement with that given by the experiment near 
the fundamental frequency (112 Hz) where the measured strains are within the sensitive 
range of the strain gages. Away from the fundamental frequency, the measured strains are 
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small compared to the sensitivity, which is approximately 1.0 /i- strain squared per Hz, and 
therefore not very accurate. As the sound pressure level is increased to 160 dB, figure 3, the 
simulation overpredicts the broadening of the spectrum and the response at high frequencies. 
It is important to restate that an equivalent linear viscous damping was used in the plate 
model, and that the nonlinear coupling to the surrounding acoustic fluid was neglected. 

• Current Model: (1) Random Excitation 

In order to obtain a pressure input spectrum similar to that of the experiments 
described above, the random data generated by the IMSL routine RNNOF is filtered using 
a quadratic narrow-band filter of bandwidth 50 to 500 Hz. Figure 4 shows that both the 
fully coupled and uncoupled models predict the same displacement frequency response of 
the plate for a low-intensity noise having a maximum sound pressure level of 120 dB. The 
spectrum obtained by the uncoupled model shows a more pronounced peak at 560 Hz, which 
corresponds to a natural frequency of the plate. The plate response is linear and is dominated 
by the first mode. Due to the high cost in CPU time of the coupled calculations, the time 
histories of the various quantities were not long enough to give smoother power spectra. 

Increasing the noise intensity to a maximum sound pressure level of 140 dB leads 
to a nonlinear plate response. The frequency content of the response shifts toward higher 
frequencies and the spectrum broadens, as shown by fig. 5. However, the shift and the 
broadening are different for the two models. The power spectrum predicted by the uncoupled 
model shows a greater shift toward high frequencies and more broadening than that predicted 
by the coupled model. The plate response at high frequencies is also higher when the 
uncoupled model is used. The difference between the coupled and uncoupled power spectra 
is similar to that observed earlier between experiments and simulation (see fig. 3). This 
result indicates that acoustic coupling at high sound pressure levels is important for the 
accurate prediction of the plate response. 

Figure 6 shows the power spectra of the pressure on both sides of the plate center. 
Notice that the pressure on the top surface is dominated by the input pressure while the 
pressure on the bottom surface is due to the vibration of the plate and has, therefore, a power 
spectrum similar to that of the response. An instantaneous pressure distribution in the Top 
domain is shown in fig. 7. Since the input is random, several peaks are observed in the pres- 
sure field superimposed on a strong low frequency. The instantaneous pressure distribution 
in the Bottom domain shows the presence of one dominant frequency corresponding to that 
of the plate response, fig. 8. 

• Current Model: (2) Harmonic Excitation 

For an excitation frequency / = 112 Hz (corresponding to the first mode of the 
plate) and an amplitude of 150 dB, the plate response is nonlinear, as shown by fig. 9. The 
power spectra of the center plate displacement show several harmonics and subharmonics 
characteristic of the nonlinear response. For this excitation level, both the fully coupled and 
uncoupled models predict the same response power spectrum (similar results axe obtained 
for lower dB levels). This result is attributed to the weak nonlinearity in the radiation field. 
Figure 10 shows the power spectra of the pressure on both sides of the plate center. The 
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power spectrum of the pressure on the top surface shows a strong peak at the fundamental 
frequency (112 Hz) and negligible harmonic content, indicative of weak nonlinearity. The 
power spectrum of the radiated pressure on the bottom surface more closely resembles to the 
response power spectrum with slightly stronger harmonics as should be expected, since the 
coupling between the plate vibration and the acoustic fluid is obtained through the vertical 
velocity. Due to transmission loss, the level of the radiated pressure on the bottom surface 
is 20 dB lower than that of the top surface. The weak nonlinearity in the pressure field 
is further evidenced by fig. 11; where the Top domain power spectra of the pressures at 
the inflow and near the center of the flexible plate are shown. The two spectra are nearly 
identical except for a negligible difference in the levels of the harmonics indicating a weak 
nonlinearity in addition to the contribution of the flexible plate. 

Increasing the level of the excitation source to 165 dB leads to a stronger contri- 
bution of the harmonics to the plate response and much less contribution from the sub- 
harmonics, fig. 12. This result is in agreement with earlier work 23 which showed that the 
window to complex dynamics is very narrow and that an increase in excitation level can 
lead to a more stable system. The two models do not predict the same displacement power 
spectrum of the plate. At high frequencies, the uncoupled model predicts a higher level of 
the response peaks. This is similar to the random excitation case, except that the peaks are 
not shifted. Figure 13 shows that the pressure fields on both sides of the plate are nonlinear. 
Both pressure power spectra show a strong harmonic content. The nonlinearity in the Top 
domain pressure field is shown in fig. 14. The levels of the harmonics on the flexible plate 
are higher due to both the plate vibrations and the nonlinear wave propagation. The latter 
point can be seen on the figure through the change in magnitude of the peak at 560 Hz 
and its harmonic 1120 Hz. On the flexible plate, the 560 Hz peak (which corresponds to a 
plate mode) is higher than its harmonic (1120 Hz). However, at the inflow the two peaks 
are nearly the same. This indicates that while the 560 Hz peak has decayed significantly its 
harmonic has decayed only slightly. This is a characteristic of nonlinear wave propagation. 
Similar results are obtained for the Bottom domain pressure field. 

5. Conclusions 

Based on the results obtained in this paper, the following conclusions can be made: 

(1) The coupling between the acoustic fluid with no meanflow and the plate vibra- 
tion is not important at low excitation levels for both random and harmonic excitations, as 
was reported by Frendi, etal. 23 

(2) The acoustic coupling is important for accurate prediction of the plate response 
at high excitation levels. The present results show that the absence of acoustic coupling in 
earlier models may be the reason for the discrepancy between numerical and experimented 
results. When acoustic coupling is accounted for, less spectral broadening and frequency 
shift is obtained. The plate response is lower at higher frequencies, which is in agreement 

n i 

with experimental observations/ 1 

(3) In the harmonic case, no shift in the peaks of the response is obtained; however, 
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the uncoupled model predicts higher levels of the higher harmonics. The results also show 
that the coupling between the acoustic fluid and the plate vibration is important only when 
both the structural response afid the acoustic radiation axe nonlinear, which is in agreement 
with previous results. 23 
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Figure 1 : Computational Domain 
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Fig. 2: Experimental and simulated power spectra of the strain for an 
overall sound pressure level of 140 dB. (Robinson, etal. **) 



Fig. 3: Experimented and simulated power spectra of the strain for an 
overall sound pressure level of 160 dB. (Robinson, etal. 
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Fig. 4: Comparison of the power spectra of the center plate displacement 
obtained by the coupled and uncoupled models for a maximum random 
excitation level of 120 dB. 
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Fig. 5: Comparison of the power spectra of the center plate displacement 
obtained by the coupled and uncoupled models for a maximum random 
excitation level of 140 dB. 
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Fig. 9: Compaxison of the power spectra of the center plate displacement 
obtained by the coupled and uncoupled models for harmonic excitation 
with a frequency of 112 Hz and a peak amplitude of 150 dB. 
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Fig. 12: Comparison of the power spectra of the center plate displace- 
ment obtained by the coupled and uncoupled models for harmonic exci- 
tation with a frequency of 112 Hz and a peak amplitude of 165 dB. 



Fig. 13: Power spectra of the surface pressure on both sides of the plate 
center. 


17 






REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No 0704-0188 


Public reDomnq Durden for this collection of information >s estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
qathering and maintaining the data needed, and completing and reviewing the collection of information Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate tor information Operations and Repons, 1215 Jefferson 
Davis Highway Suite 1204 Arlington. VA 22202-4302. and to the Office of Management and Budget. Paperwork Reduction Project (0704-0188), Washington, DC 20503 * 


1. AGENCY USE ONLY (Leave blank ) 


2. REPORT DATE 
March 1993 


3. REPORT TYPE AND DATES COVERED 
Technical Memorandum 


4. TITLE AND SUBTITLE 

On the Effect of Acoustic Coupling on Random and Harmonic Plate Vibrations 

5. FUNDING NUMBERS 
WU 505-63-50-10 

6. AUTHOR(S) 

A. Frendi 
J. H. Robinson 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 
NASA Langley Research Center 
Hampton, VA 23681-0001 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

9. SPONSORING /MONITORING AGENCY NAME(S) AND ADDRESS(ES) 
National Aeronautics and Space Administration 
Washington, DC 20546 

10. SPONSORING /MONITORING 
AGENCY REPORT NUMBER 

NASA TM 107740 

11. SUPPLEMENTARY NOTES 

A. Frendi: Analytical Services & Materials, Inc.; Hampton, VA 
Robinson: NASA Langley Research Center, Hampton, VA 

12a. DISTRIBUTION /AVAILABILITY STATEMENT 

Unclassified — Unlimited 
Subject Category— 71 

12b. DISTRIBUTION CODE 


13. abstract (Maximum 200 words) gffeg 0 f acoustic coupling on random and harmonic plate vibrations is studied using 

two numerical models. In the coupled model, the plate response is obtained by integration of the nonlinear plate equation 
coupled with the nonlinear Euler equations for the surrounding acoustic fluid. In the uncoupled model, the nonlinear plate 
equation with an equivalent linear viscous damping term is integrated to obtain the response of the plate subject to the same 
excitation field. For a low-level, narrow-band excitation, the two models predict the same plate response spectra. As the 
excitation level is increased, the response power spectrum predicted by the uncoupled model becomes broader and more 
shifted towards the high frequencies than that obtained by the coupled model. In addition, the difference in response 
between the coupled and uncoupled models at high frequencies becomes larger. When a high intensity harmonic excitation 
is used, causing a nonlinear plate response, both models predict the same frequency content of the response. However, the 
level of the harmonics and subharmonics are higher for the uncoupled model. Comparisons to earlier experimental and 
numerical results show that acoustic coupling has a significant effect on the plate response at high excitation levels. Its 
absence in previous models may explain the discrepancy between predicted and measured responses. 


15. NUMBER OF PAGES 

19 


14. SUBJECT TERMS 

Fluid-Structure Coupling, Acoustic Damping, Acoustic Radiation, Random Vibration, Vibration, 

Nonlinear, Plates 1 16 - PR,CE C0DE 

A03 


17. SECURITY CLASSIFICATION I 18 SECURITY CLASSIFICATION 119. SECURITY CLASSIFICATION I 20. LIMITATION OF ABSTRACT 

OF REPORT I OF THIS PAGE I OF ABSTRACT 


Unclassified 


NSN 7540-01-280-5500 


ndassified 




Standard Form 298 (Rev 2-89) 

Pn?^cnb ,a d bv ANSI Std 















